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FOR ATTITUDE DETERMINATION USING 
GLOBAL POSITIONING SYSTEM SIGNALS 

John L. Crassidis,* F. Landis Markley,* E. Glenn Lightsey* 

In this paper, a new motion-based algorithm for GPS integer 
ambiguity resolution is derived. The first step of this algorithm 
converts the reference sightline vectors into body frame vectors. 

This is accomplished by an optimal vectorized transformation of 
the phase difference measurements. The result of this 
transformation leads to the conversion of the integer ambiguities to 
vectorized biases. This essentially converts the problem to the 
familiar magnetometer-bias determination problem, for which an 
optimal and efficient solution exists. Also, the formulation in this 
paper is re-derived to provide a sequential estimate, so that a 
suitable stopping condition can be found during the vehicle 
motion. The advantages of the new algorithm include: it does not 
require an a-priori estimate of the vehicle’s attitude; it provides an 
inherent integrity check using a covariance-type expression; and it 
can sequentially estimate the ambiguities during the vehicle 
motion. The only disadvantage of the new algorithm is that it 
requires at least three non-coplanar baselines. The performance of 
the new algorithm is tested on a dynamic hardware simulator. 

INTRODUCTION 

The utilization of phase difference measurements from Global Positioning System 
(GPS) receivers provides a novel approach for three-axis attitude determination and/or 
estimation. These measurements have been successfully used to determine the attitude of 
air-based , 1 space-based , 2 ' 3 and sea-based 4 vehicles. Since phase differences are used, the 
correct number of integer wavelengths between a given pair of antennas must be found. 
The determination of the integer ambiguities can either be accomplished by using “static” 
(motionless) or “dynamic” (motion-based) techniques. The ambiguities essentially act as 
integer biases to the phase difference measurements. Once the integer ambiguities are 
resolved, then the attitude determination problem can be solved . 5 

The static method finds a solution that minimizes the error residual at a specific time 
by searching through an exhaustive list of all possible integers and rejecting classes of 
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solutions when the residual becomes too large. 6 Refinements can be made to the solution 
by restricting the search space with knowledge of a-priori information, such as the 
maximum tilt the baseline should encounter. 7 Static methods generally rely on solving a 
set of Diophantine equations. 8 The appeal of these methods is that they provide an 
“instantaneous” attitude solution, limited only by computation time, and are well suited 
for the short baselines. However, the minimum residual does not guarantee a correct 
solution in the presence of noise. 1 In fact, it is possible that static methods can report a 
wrong solution as valid, especially when some of the calibration information, such as line 
bias, is incorrect. This lack of integrity can cause significant problems if the sensor 
output is used to control a high bandwidth actuator, such as gas jets on a spacecraft. 
Another consideration is that static methods sometime require that the antenna array must 
be within a defined angle (typically 30 degrees) of a reference attitude, which is often true 
for ground-based applications, but is less likely for space-based applications. Also, 
structural flexibility in the baselines may lead to erroneous solutions. All of the 
aforementioned limitations imply that static methods, while attractive because of their 
fast solutions, are not totally acceptable for general purpose applications. 

The other technique for resolving integer ambiguities involves collecting data for a 
given period of time and performing a batch solution, in which the integer terms remain 
constant over the collection period. This technique relies on the fact that a certain amount 
of motion has occurred during the data collection, either from vehicle body rotation or 
GPS line of sight motion. The main disadvantage of this technique, compared to static 
approaches, is that it takes time for the motion to occur, which may be on the order of 
several minutes. Another consideration is that a potentially significant amount of 
memory is required for the storage of the batch data collection. But, motion-based 
techniques also have significant advantages over static methods. Most importantly, 
motion-based techniques are inherently high integrity methods because there are 
numerous checks that can be implemented into the solution before it is accepted. These 
include using statistical checks applied to error residuals, matrix condition number 
checks, and using the closeness of the computed floating-point “integers” to actual 
integers as a check. The probability of an erroneous solution being reported as valid can 
be made as small as desired by appropriately setting the thresholds on these integrity 
checks. For these reasons, motion-based techniques have been more widely used for on- 
board applications. 

Traditional motion-based techniques of integer ambiguity resolution rely on the fact 
that either GPS line of sight motion or vehicle motion dominates the changes in 
differential carrier phase measurements. Cohen 9 developed an algorithm, known as 
“quasi-static” integer resolution, that can be used when the GPS line of sight motion and 
the vehicle rotation both account approximately evenly for the differential carrier phase 
measurement changes. This algorithm can be adapted to almost any vehicle motion, slow 
or fast, simply by varying the sample rate and the data collection time. The quasi-static 
method solves a collection of differential phase measurements for a single attitude 
estimate and then considers perturbations to the initial estimate at each measurement 
epoch to produce a time varying batch solution to the data. Although this is a widely used 
algorithm, there are certain disadvantages. First, an a-priori attitude estimate must be 
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given. Second, the algorithm is an iterative batch estimator that may produce erroneous 
estimates, depending on the accuracy of the a-priori attitude estimate. Finally, if a large 
number of samples in the data collection are required to observe the motion, large-order 
matrix inversions may be required. Another method (Ref. 10) performs a minimization 
on three Euler-angle attitude parameters independent of each other, followed by 
determining the integers. This approach has been shown to provide better convergence 
than Cohen’s method and works well for non-coplanar baselines; however, singular 
conditions can exist at various attitude rotations and a significant amount of vehicle 
motion may be necessary for a solution. 

In this paper, a new motion-based algorithm is derived. The main advantages of the 
new algorithm over the prior methods include: (i) it resolves the integer ambiguities 
without any a-priori attitude knowledge, (ii) it requires less computational effort, since 
large matrix inverses are not needed, and (iii) it is non-iterative. The only disadvantage 
of the new algorithm is that it requires at least three non-coplanar baselines. The 
algorithm is first shown as a batch solution, and then shown as a sequential solution. A 
covariance expression is also derived which can be used to bound the integer solution so 
that a sufficient integrity check for convergence can be developed. This is extremely 
useful in the sequential formulation, since the solution can be found as the motion occurs, 
rather than taking a batch solution at a specific data collection interval. For these reasons, 
the new algorithm provides an attractive method for real-time ambiguity resolution. 

The organization of this paper proceeds as follows. First, the concept of the GPS 
phase difference measurement is introduced. Then, a brief review of Cohen’s quasi-static 
method is shown, and limitations and computational aspects of this algorithm are 
discussed. Next, the new motion-based algorithm is derived. The conversion of the GPS 
sightline vector into the body frame is first reviewed. Then, the batch solution used to 
resolve the integer ambiguities is derived, followed by the sequential solution. Finally, 
the new algorithm is validated by using an actual GPS receiver with a hybrid dynamic 
simulator to simulate the vehicle motions of a low-altitude Earth-orbiting spacecraft. 

GPS SENSOR MODEL 

In this section, a brief background of the GPS phase difference measurement is 
shown. The main measurement used for attitude determination is the phase difference of 
the GPS signal received from two antennas separated by a baseline. The wavefront angle 
and wavelength are used to develop a phase difference, as shown in Figure 1. The phase 
difference measurement is obtained by 

bi cos 6 = A( A<f> - n) (1) 

where b[ is the baseline length (in cm), 0 is the angle between the baseline and the line of 
sight to the GPS spacecraft, n is the number of integer wavelengths between two 
receivers, A<f> is the phase difference (in cycles), and X is the wavelength (in cm) of the 
GPS signal. The two GPS frequency carriers are LI at 1575.42 MHz and L2 at 1227.6 
MHz. As of this writing, non-military applications generally use the LI frequency. The 
phase difference can be expressed by 

A<f> = b T As+n (2) 
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where seR is the normalized line of sight vector to the GPS spacecraft in a reference 
frame, beR is the baseline vector (in wavelengths), which is the relative position vector 
from one receiver to another, and A e /? 3x3 is the attitude matrix, which is an orthogonal 
matrix with determinant 1 (i.e., A T A = I 3*3). The measurement model is given by 

M>ij = kj A sj + rijj + wy ( 3 ) 

where A 0 y denotes the phase difference measurement for the i* baseline and f' sightline, 
and represents a zero-mean Gaussian measurement error with standard deviation 
which is 0.5 cm/ A = 0.026 wavelengths for typical phase noise. 9 


To GPS 



Figure 1 GPS Wavelength and Wavefront Angle 

QUASI-STATIC APPROACH 

Cohen’s quasi-static method 9 is a motion-based technique that begins by t akin g 
measurements for k = 1 to L (“measurement epochs”) to which a single attitude solution 
will be determined. At each epoch it is assumed that M baselines exist and N 
sightlines. The measurement model is linearized by assuming a small perturbation about 
a nominal attitude Aq and an assumed set of integer phases (no) - , so that 

A = 4 (/ 3x3 +[j»x]) ( 4 ) 

where 73x3 is a 3 x 3 identity matrix, 56 is assumed to be a small angle rotation, and 
[<5ff x] is a cross product matrix with 
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(5) 


0 -fl3 <*2 

03 0 -a\ 

-02 <2] o 

Substituting Equation (4) into Equation (3) for all available measurements yields 

r iMfex] 

S# = [ A*-A*J = : / 

-[* ' 

where 7 is a quasi-identity matrix with possible zeros along the diagonal where states 
have been removed at various measurement epochs, and 

’ Ml 1 (*-«o) n 

: , &= \ , (A0q)„ = bf Aqs j + (n 0 ) z y (7) 

Equation (6) is a set of MN equations for (3 + MN) states. Allowing perturbations at all 
epochs leads to 

3p(l)l \H(l) 0 ••• 0 7(1) lr 30(1) 

30(2) _ 0 H{ 2) ! 7(2) i 

: “ i \ 0 : S0(L) 

S0(L)\ 0 ••• 0 H(L) 7(L)jL 

This compact representation has LMN rows and (3Z, + MN) states. In principle, the 
integers and the attitude of the vehicle of the vehicle at each measurement epoch may be 
found by applying an iterative linearized least-squares approach using Equation (8) and 
updating the nominal attitude using Equation (4). 

The quasi-static method has been successfully implemented to resolve the integer 
ambiguities on an actual system (Ref. 1), and works extremely well when a fairly accurate 
a-priori attitude is known, and significant vehicle motion is present. However, this 
approach has a number of disadvantages. First, if the a-priori attitude estimate is poorly 
known, then the solution may never converge (even if the integers are known exactly). 
Second, not only does this algorithm require a good a-priori guess, but requires fairly 
accurate attitude estimates at all measurement times. The reason for this is that as time 
increases, the perturbations to the a-priori attitude guess may become too large for the 
solution to converge. This may be overcome by augmenting the state equations to 
include a constant, but unknown, body rate that is also estimated. Finally, a (3 L + MN) 
matrix inverse is required, which may cause computational problems. Many of these 
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problems may be overcome by developing an algorithm that is independent of any 
attitude information. 


NEW ALGORITHM 

In this section a new algorithm to resolve the integer ambiguities is shown. The main 
advantage of this algorithm is that it is attitude independent. First, a conversion of the 
sightline vectors into the body frame is shown. This converts the problem into the 
familiar magnetometer-bias problem. Then, a batch solution for this problem is shown, 
followed by a sequential approach. 

The new algorithm begins by determining the sightline vector in the body frame, 
denoted by s = As. This is accomplished by minimizing the following loss function 1 1 



If at least three non-coplanar baselines exist, the minimization of Equation (9) is 
straightforward and leads to 



-l 


Cj=Bj 


M 

,=i w u 


M 

/=i m ii 


(10a) 

(10b) 

(10c) 

(lOd) 


The computed sightline in the body frame is related to the sightline vector in the reference 
frame by 

Sj=ASj+Cj+£j ( 11 ) 


where Cj is a constant bias since the baselines are assumed constant, and e j is a zero- 

mean Gaussian process with covariance Rj = BJ l . Again, the inverse in Equation (10) 

exists only if three non-coplanar baseline vectors exist. 

The next step is to use an attitude-independent method to find the phase-bias vector 
cj. Doing this for each sightline gives us all the sightlines in both the body frame and the 

reference frame. The explicit integer phases are not needed for this solution, but it is 
important to check that they are close to integer values, as mentioned in the Introduction. 
In the general case, the explicit integer phases can be found from the attitude solution. 
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The three-baseline case (M = 3) is simpler, for in this case Equation (10c) can be inverted 
to give 


n ij —i -j 


( 12 ) 


With more than three baselines, however, Equation (10c) does not have a unique solution 
for Cj, so the M integer phases for sightline Sj cannot be found from Cj alone. We will 

consider the three-baseline case, which is the most common in practice. If more baselines 
are available, we are always free to select a three-baseline subset. Then, after the integer 
phases have been determined, a refined attitude estimate can be computed using all 
baselines (i.e., three baselines are sufficient to determine an attitude, which may then be 
used to resolve the integers corresponding to the other baselines). 

To eliminate the dependence on the attitude, the square of Equation (1 1) is computed, 
so that 





(13) 


Next, the following effective measurement and noise are defined 



(14a) 

(14b) 


Then, the effective measurement can be written as 


Zj — 2 5 j * Cj ||c 7 - 1 + V j 


( 15 ) 


Alonso and Shuster (Ref. 12) showed that vj is approximately Gaussian for small Sj with 
mean given by 

jjj = £{vjj = -trace|/?j} (16) 


and variance given by 

aj =E{v]}-M 2 j = 4 (sj-CjfRjlij-Cj) 

The negative-log-likelihood function for the bias is given by 

= z #)- 2 M*) £ HM| 2_ ^(*>] +1 °g o- 2 j(k) + \og2x 


(17) 


(18) 


The symbol k denotes the variable at time t k . The maximum-likelihood estimate for Cj, 
denoted by c*-, minimizes the negative-log-likelihood function, and satisfies 
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(19) 



The minimization of Equation (18) is not straightforward since the likelihood function 
is quartic in Cj. A number of algorithms have been proposed for estimating the bias (see 

Ref. 12 for a survey). The simplest solution is obtained by scoring, which involves a 
Newton-Raphson iterative approach. Another approach avoids the minimization of a 
quartic loss function by using a “centered” estimate. A statistically correct centered 
estimate is also derived in Ref. 12. Furthermore, Alonso and Shuster show a complete 
solution of the statistically correct centered estimate that determines the exact marimnm 
likelihood estimate c*j. This involves using the statistically correct centered estimate as 

an initial estimate, and iterating on a correction term using a Gauss-Newton method. 
Although this extension to the statistically correct centered estimate can provide some 
improvements, this part is not deemed necessary for the GPS problem since the estimated 
quantity for is rounded to the nearest integer. 

Batch Solution 


In this section the statistically correct centered estimate algorithm (see Ref. 12 for 
details) and its application to the integer ambiguity problem are shown. First, the 
following weighted averages are defined 


1 L \ 


k=l 

L 


= v#) ’ ^ 


( 20 ) 


where 


1 



z 


Jt=l 



( 21 ) 


Next, the following variables are defined 

Zj(k) = Zj(k)-Zj, lj(k) = Sj(k)-lj, Vj(k) = Vj(k)-vj, Jij{k) = jUj(k)-Jij (22) 


The statistically correct centered estimate now minimiz es the following loss function 



-^[zj(k)-2lj(k)-Cj- M j(k)] 2 


(23) 


which is now a quadratic function in Cj. The minimiza tion leads directly to 
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£ 'i = p i < 24 > 

where the estimate error covariance is given by 

' L T 1 

P J= Z^27 77 4 M*)lJ(*) ( 25 ) 

_k=l a J {K) 

The ambiguity for the I th baseline and 7 th sightline can be resolved by rounding the 
following to the nearest integer 

n ij - kf £*j (26) 

The integer error covariance, denoted by Qj, can be shown to be given by 

Qij -tfPjki (27) 

Equation (27) can be used to develop an integrity check for the algorithm. For example, a 
suitable criterion can be developed from a three-sigma bound using 3 yfQjj- 

Sequential Formulation 

This section expands upon the batch solution so that a sequential estimate of the 
integers can be found. The main advantage of a sequential formulation is that the 
convergence (integrity) check can be made on-the-fly (i.e., in real-time). The covariance 
in Equation (25) to be expanded to the L + 1 time point, so that 

L 

P7\l + 1) = y -4— 4 s , (k)s? (Jfc) + 4 s , (L + 1) s] (L + 1) 

rM(*) 7 o 2 al+ 1) ~ J J 

*=1 JK ’ ]K ’ (28) 

= P" 1 (L) + -^ 4?,(L + l)r[(L + l) 

J aj(L + 1 ) ~ J J 

From the matrix inversion lemma, 13 the following sequential formation for the covariance 
is developed 

Pj(k + l) = Kj(k)Pj(k) (29) 

where 

Kj(k) = I - Pj(k)lj(k + 1) ?J (k + 1) Pj(k)lj(k + 1) + ±oj(k + 1)1 Jj (k + 1) (30) 

In order to derive sequential formulas for the quantities in Equation (20), first consider 
the following identity 

273m = i?(Z) ^ = fj<L> (31) 
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Expanding out this expression using L - 1 points in the summation yields 


&A L ~ >) Zy(i ^(I-l) + ^(I) 


and so 


?/(£) = 


gj(L); / (L-l)+gj(Z.-l);,(L) 


aj(L)+oj(L-l) 


(34c) 


Therefore, the following sequential expressions for the quantities in Equation (20) are 
given 

+ 1 )= m &j(k+ + + + ^ (34a) 

s_j (£ + 1) = + ^ + (£ + l)ly(fc)+0y (*)!/(* + l)] (34b) 

+ 1 ) = q .2(^ + 1) 1 + -2 w H + 1 W*) + WVjtt + 1)] (34c) 

where 

ct^(A: + 1) cry(fc) cPj(k + \) ^ ^ 

The estimated bias in Equation (24) can also be found in a similar mann er, so that 

c*{k + 1) = Kj{k)c*{k) + ~ [zj(k + 1) - Jij{k + 1)]2 Pj{k + 1 )lj(k + 1) (36) 

4- 1 j 

Since the baselines are constant, Equations (26) and (27) can be used directly to 
determine the sequential integer value and error covariance, given by 


where 


n ij(k) - kf c*j(k) 
Q i j(k) = bfp j (k)b i 


(37a) 

(37b) 


The complete solution proceeds as follows. First, use Equations (10b) and (lOd) to 
convert the sightline vectors into the body frame. Then, perform an initial batch solution 
using Equations (20)-(25) in order to initialize the sequential routine (an accurate initial 
estimate is not required as will be seen in the results section). Then, perform a sequential 
estimate for the integers using Equations (29), (30), and (34)-(37). Finally, continue until 
the covariance in Equation (37b) is below a pre-specified value. 

There are many advantages of the new algorithm. First, the algorithm is fully 
autonomous (i.e., it requires no a-priori information such as an a-priori attitude guess). 
Second, the largest matrix inverse is of a 3x3 matrix, which makes the algorithm 
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computationally efficient and stable. Third, it is non-iterative, which makes it suitable for 
a sequential formulation. This has a significant advantage since the convergence can be 
checked during the actual motion in the vehicle. Finally, the integers for other sightlines 
can be easily resolved by calling the same subroutine. Therefore, the algorithm can easily 
be implemented using all available sightiines, and attitude determination can begin once 
the integers corresponding to two sightlines have been resolved. For these reasons, the 
new algorithm provides an attractive approach to resolve the integers. 

HARDWARE SIMULATION AND RESULTS 

A hardware simulation of a typical spacecraft attitude determination application was 
undertaken to demonstrate the performance of the new algorithm. For this simulation, a 
Northern Telecom 40 channel, 4 RF output STR 2760 unit was used to generate the GPS 
signals that would be received at a user specified location and velocity. The signals are 
then provided directly (i.e., they are not actually radiated) to a GPS receiver that has been 
equipped with software tracking algorithms that allow it operate in space (see Figure 2). 

Inputs Visible Satellites 

GPS Constellation Doppler Shift 



Performance Position 

Velocity 

Attitude 


Figure 2 Hardware Simulation Block Diagram 

The receiver that was used was a Trimble TANS Vector; which is a 6 channel, 4 RF 
input multiplexing receiver that performs 3-axis attitude determination using GPS earner 
phase and line of sight measurements. This receiver software was modified at Stanford 
University and NASA-Goddard to allow it to operate in space. This receiver model has 
been flown and operated successfully on several spacecraft, including: REX-II, OAST- 
Flyer, GANE, Orbcomm, Microlab, and others. 
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The simulated motion profile was for an actual spacecraft, the Small Satellite 
Technology Initiative (SSTI) Lewis satellite, which carried an experiment to assess the 
performance of GPS attitude determination on-orbit. Although the spacecraft was lost 
due to a malfunction not related to the GPS experiment shortly after launch, this motion 
profile is nonetheless very representative of the types of attitude determination 
applications. The orbit parameters and pointing profile used for the simulation are given 
in Table 1. 


Table 1 SSTI Lewis Orbit parameters 


Semimajor axis (a) 

6901.137 km 

Inclination (i) 

97.45 deg 

Right Ascension of Ascending Node (RAAN) 

-157.1 deg 

Eccentricity (e) 

0.0001 

Pointing profile 

Earth pointed 

Launch date 

August 22, 1997 


The simulated SSTI Lewis spacecraft has four GPS antennas that form three 
baselines. The antenna separation distances are 0.61 m, 1.12 m, and 1.07 m, respectively. 
One antenna (in baseline 3) is located 0.23 m out of plane (below) the other three 
antennas. On the spacecraft, the antennas are mounted on pedestals with ground planes to 
minimize signal reflections and multipath. For the simulation, the signal was provided to 
the GPS receiver without multipath noise. The baseline vectors in wavelengths are given 
by 


‘ 2.75 ' 


"0.00' 


-3.93" 

1.64 

K> 

II 

6.28 

II 

CO 

^1 

3.93 

-0.12 


-0.17 


-1.23 


Line biases are first determined before the new algorithm is tested to resolve the integer 
ambiguities. The GPS raw measurements are processed at 1 Hz over a forty minute 
simulation. During the simulated run, a minimum of three visible GPS are given at ail 
times. Also, there are a number of eight minute spans when two of the same (in tim e) 
sightlines are available for the ambiguity resolution algorithm. Again, in practice, all 
available sightlines should be processed simultaneously, but with three baseline vectors 
only two simultaneously available sightlines are required to determine the attitude of the 
vehicle. 

As mentioned previously, the first step in the algorithm involves using the baselines 
and phase difference measurements to convert the sightline vector into the body-frame, 
using Equations (10b) and (lOd). Then, a small batch run is used to initialize the 
sequential routine. For this case, only 5 seconds of data was required to perform the 
initialization. Again, only two sightlines are required to determine the attitude. 
Sequential error results (i.e., actual integer minus the computed values without rounding) 
for the first sightline are shown in Figure 3. The integer error can be found be rounding 
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the values in Figure 3 to the nearest integer. Clearly, for this case, the integer ambiguities 
have been resolved even before the sequential process begins (i.e., within 5 seconds). A 
plot of the 'i^Qy values is shown in Figure 4 for the first sightline (a suitable integrity 

check is given when 3^fQ~ is below 0.5). Clearly, the integrity check shows that the 

ambiguities are resolved within 5 minutes. Note, that this is a sufficiency test (i.e., the 
integers may be resolved well before 5 minutes, which is seen in this case). A plot of the 
errors for the second sightline is shown in Figure 5. For this case, all of the ambiguities 
have been resolved within 30 seconds (the error value corresponding to the second 
baseline goes below -0.5 before 30 seconds). A plot of integrity check for the second 
sightline is shown in Figure 6. The integrity check shows that the ambiguities are 
resolved within 7 minutes. This hardware simulation of a spacecraft clearly demonstrates 
that the new algorithm presented in this paper provides an accurate method to resolve the 
integer ambiguities with even slight vehicle motion. 




Figure 3 Errors for First Sightline Figure 4 Integrity Check for First Sightline 



Figure 5 Errors for Second Sightline 



Figure 6 Integrity Check for Second Sightline 
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CONCLUSIONS 

In this paper, a new algorithm was developed for GPS integer ambiguity resolution. 
The new algorithm has several advantages over previously existing algorithms. First, the 
algorithm is attitude independent so that no a-priori attitude estimate (or assumed vehicle 
motion) is required. Second, the algorithm is sequential so that it may be implemented in 
real-time. Also, a suitable integrity check can be used to determine when the determined 
values have converged to the correct values. Finally, the algorithm is computationally 
efficient since only a 3x3 matrix inverse is required, and the same subroutine can be 
used on different sightlines. The only disadvantage of the new algori thm is that it 
requires at least three non-coplanar baselines. The algorithm was tested using a GPS 
hardware simulator to simulate the motions of a typical low-altitude Earth-orbiting 
spacecraft. Results indicated that the new algorithm provides a viable and attractive 
means to effectively resolve the integer ambiguities. 
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